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ABSTRACT 

We describe a method from which cosmology may be constrained from the 2QZ 
Survey. By comparing clustering properties parallel and perpendicular to the line of 
sight and by modeling the effects of redshift space distortions, we are able to study 
geometric distortions in the clustering pattern which occur if a wrong cosmology is 
assumed when translating redshifts into comoving distances. Alcock and Paczyhski 
pointed out that this technique was particularly sensitive as a test for a large cosmo- 
logical constant, A. 

Using mock 2QZ catalogues, drawn from the Hubble Volume simulation, we find, 
that there is a degeneracy between the geometric distortions and the redshift-space 
distortions, parameterised by /3qso(z), that makes it difficult to obtain an unambigu- 
ous estimate of f2 m (0), the matter density parameter, from the geometric tests alone. 
This is in agreement with the conclusions of Ballinger et al. However, we demonstrate 
a new method to determine the cosmology which works by combining the above geo- 
metric test with a test based on the evolution of the QSO clustering amplitude, which 
has a different dependence on (3qso(z) and fi m (0). In the analysis of the Hubble Vol- 
ume mock catalogues we find that we are able to break the degeneracy between r2 m (0) 
and (3qso(z) and that independent constraints to ±20% (ltr) accuracy on r2 m (0) and 
±10% (la) accuracy on (3qso(z) should be possible in the full 2QZ survey. 

Finally we apply the method to the 10k catalogue of 2QZ QSOs. The smaller 
number of QSOs and the current status of the Survey mean that a strong result on 
cosmology is not possible but we do constrain (3qso{z) to 0.35 ± 0.2. By combining 
this constraint with the further constraint available from the amplitude of QSO clus- 
tering, we find tentative evidence favouring a model with non-zero Oa (0), although an 
^m(0) = 1 model provides only a marginally less good fit. A model with Qa(0) = 1 
is ruled out. The results arc in agreement with those found by Outram et al. using a 
similar analysis in Fourier space. 
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1 INTRODUCTION 

One of the challenges that still faces cosmologists is to de- 
termine the values of the underlying cosmological parame- 
ters of the Universe, such as the density parameter, fi m (0), 
and the energy density associated with the cosmological 
constant, ^a(0). During the 1990's, it became clear that 



there was a problem with the favoured Einstein-de Sitter, 
n m (0) = l model. Hubble's constant was being measured as 
Ho ~ 70 km s _1 Mpc _1 which gives an age for the Uni- 
verse of ~ 10 Gyr if fi m (0) = l. Globular clusters were es- 
t imated to be older than the Universe itself at 16±2 Gyr 
( Etenzini et al.| 1996), although recent estimates of the ages 
of Globular clusters are slightly lower, e.g. 11.8±2.1 Gyr 
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( Gratton ct al.||l997 ). However, inflationary theory, in its 
simplest form, implies a flat universe. To comply with this 
requirement, a cosmology where f2 m (0) + f2A(0) = l was sug- 
gested as an alternative to the Einstein-de Sitter model 
(Peebles 1984). Invoking a cosmological constant, A term 
(£2a(0) = Ac /(3-STo)), provides a solution to the problem 
of the age of the Universe, since with lower values of fi m (0) 
(and fi m (0) + fiA(0) = l), the age of the Universe is increased. 
A cosmology with non-zero A also matches the amplitude of 
fluctua tions on both COBE ( gmoot et al. 1992 ) and cluster 
scales ( Eke, Cole fc Frcnk 1996 ). However, the physics un- 
derpinning such a small A term remains very unclear (e.g. 
Witten 2000). 

One of the most elegant ways of testing for the exis- 
tence of a non-z ero A term was put forward by Alcock & 
Paczyngku (1979). They suggested that if the clustering of 
galaxies or QSOs, parallel and perpendicular to the line of 
sight, was computed assuming an Einstein-de Sitter cosmol- 
ogy but the Universe had a different cosmology, for example 
one which had Oa(0) 7^ 0, then a squashing or elongation 
of the clustering pattern would occur more in one direction 
than in the other. This is because assuming the wrong cos- 
mology affects the clustering pattern parallel to the line of 
sight differently to that perpendicular to the line of sight. 
The cosmology assumed for measurements of the clustering 
can have any value of f2 m (0) and £2a(0) but if these values 
differ from the underlying cosmology of the Universe, distor- 
tions are introduced into the clustering pattern. The only 
assumption required is that the clustering in real space is on 
average spherically symmetric. 

Clustering statistics measured from galaxy or QSO red- 
shift surveys are, however, not measured in real space and 
so they have redshift space distortions imprinted on them. 
On small scales, virialized clusters appear elongated along 
the line of sight in redshift space. These are the so called 
'fingers of God'. They arise because the redshift due to cos- 
mological expansion and redshift due to peculiar motions 
cannot be separated. The redshift is converted into a dis- 
tance assuming that there are no peculiar motions, which 
distorts the shape of the cluster. On large scales, coherent 
infall squashes over-densities along the line of sight in red- 
shift space. This causes a boost in the amplitude of the red- 
shift space correlation function on linear scales, as compared 
to the real space correlation function. Th is is characterised 
by the parameter j3 = f2^ 6 /6 ( Kaiser 1987 ). In order to mea- 



sure properly any effect of A, these two contributions have 
to be accurately accounted for in any model of the cluster- 
ing. Models for the redshift spac e distortions have appeared 
many t imes in the literature (e.g. Bcaiictaljll983|, baningcr 
Peacocl : fe Heavens 199( , Matsubara fe Suto|[199G| , [Popowsk 
et al. 1 )9£, Ratcliffe et al. 1998c), and we describe some of 
these models in section 14 

Either galaxy or QSO surveys can, in theory, be used 
to detect the distortion introduced into the clustering pat- 
tern through incorrect assumptions for the cosmology. How- 
ever, wide angle galaxy surveys, including the 2dF Galaxy 
Redshift Survey, essentially probe the Universe at z = 0, 
where assuming the incorrect cosmology would make lit- 
tle difference to the measured clustering pattern. Cluster- 
i ng from pencil beam surveys, such as the CNOC survey 
( Yee et al. 200Q ), cannot be measured over a wide range of 
scales. QSO surveys sample space more sparsely than cur- 



rent galaxy surveys but probe clustering out to high redshifts 
over a wide area, potentially allowing the distortions in the 
clustering to be detected. 

The aims of this paper are two fold. First we use the 
mock catalogues to test how well f2 m (0), £2a(0) and Pqso(z) 
will be constrained from the 2dF QSO Redshift Survey 
(2QZ). At this stage, we only consider flat cosmologies with 
fl m (0) + flA(0) = l, as expected from current CMB experi- 



ments (Balbi et al. 2000, de Bernardis et al. 2000) and in- 



flationary theory. The method can, however, be extended to 
allow for open or closed cosmologies. We also test whether 
tighter constraints on the cosmology can be found if the 
results from £(<r, n) are combined with results from the evo- 
lution of the QSO-mass bias . Second, we apply t he method 
to the 2QZ 10K Catalogue ( Croom et al, 2001b ) to obtain 
preliminary results on cosmology and demonstrate the fea- 
sibility of this method when applied to real data. 

The outline of this paper is as follows. In section fel we 
describe both the 2QZ Survey and the Hubble Volume simu- 
lation, from which the mock catalogues are constructed. In 
section H, we discuss how the correlation function parallel 
and perpendicular to the line of sight, f(<7, 71"), is measured 
from either the simulation or the 2QZ Survey and in sec- 
tion W we discuss how we model £(<r, n). In section feL we 
outline our method for detecting cosmology and in section 
r& we present our results. In section M, we consider whether 
other constraints can be usefully combined with the results 
from £(g, 7 r). In Section p we a pply the analysis to the 10k 
catalogue (Croom et al. pOOlb) and finally in section H we 



draw our conclusions. 



2 2QZ SURVEY AND MOCK CATALOGUES 

2.1 2QZ 

The 2QZ redshift survey aims to measure the redshifts of 
25,000 QSOs over a redshift range of 0.3< z < 3. The spec- 
tra of the QSOs will be measured using the 2dF instrument 
on the Anglo- Australian Telescope (AAT). This instrument 
allows up to 400 spectra to be obtained simultaneously, al- 
lowing a large number of QSOs to be observed in a relatively 
short period of time. QSOs are selected in a homogeneous 
manner via ubjr multi-colour selection and the survey wi ll 
be more than 90% complete to z ~2.2 ( [Boyle et ah| ^000| ). 
The spectra of the objects are obtained using the 2dF in- 
strument on the AAT. T he data is reduced using the 2dF 
pipeline reduction system (Bailey and Glazebrook 1999) and 



objects are identi fied as QSOs by an automated procedure 
known as AUT0Z ( |Miller et al] fcOO l[ ). AUT0Z also identifies 
the QSO redshifts; these are then visually confirmed by two 
independent observers. 

Two strips of the sky are being observed, each mea- 
suring 75° x 5°. One is centred at S = -30°, with 21 h 40 m 
<C a <. to 03 h 15 m , close to the Southern Galactic Pole, while 
the other is centred at 8 = 0° with 09 h 50 m < a < to 14 h 50 m 
in the North Galactic Cap. The total area will be around 740 
square degrees. 

For more details see 
aT 



ioom et al 



(200C1 



Croom et al 



Cro om et al 

( EooibQ , 



(|2001aD, |Hoyle et al 



(|!998d), Boyle et 



Smith et al 



(2001) 



( pooiD , 

and also 



www . 2df quasar . org. 
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2.2 The Hubble Volume Simulations 

The Hubble Volume Simulations r] are a series of simulations 
run by the Virgo Consortiu m and to date are th e larg est N- 



body s imula tions ever (see Jenkins et al. (1998) and Frcnk 



et al, (2000) for more details on Virgo consortium simula- 



tions). In order to make realistic mock catalogues, a simu- 
lation where the dark matter particles are output along an 
observer's past lightcone is required. The ACDM simulation 
has one deep lightcone output which covers an area of 15 x 
75 degrees and extends out to z ~ 4. 

The parameters of the simulation are f2 m (0)=0.3, 
S1a(0)=0.7, H o =70 km s Mpc -1 and the normalisation, 
o~$, is 0. 9, consistent with the abundance of hot X-ray 
clusters (White, Efstathiou fc Frcnk 1993) and with the 



level of anis otropics in the CMBR found by COBE (Smoot 
et al. 1992). The input power spectrum was calculated 
using CMBFAST flSeljak fc Zaldarriaga| |l996| ) with the 



parameters listed above but assuming that f2 b ( O )=0.04 
and ficDM(o) = 



b(0) = 

:0.26. This changes the shape of the in- 



O. m (0) h exp(-r2 b(0 )[v2ft + 
as compared to T = 



199! 



put power spectrum to I 
On(0)]/fi m (0))=0.17 (|Sugiyam; i 
fi m (0)/i = 0.21. 

One billion mass particles are contained within a cube 
that is 3,000/i~ Mpc on a side. One of the vertices was 
chosen to be the observer and the long axis of the light- 
cone was oriented along the maximal diagonal. The light- 
cone therefore extends to a depth of ~5,000/i -1 Mpc, which 
corresponds to z ~4 in the ACDM cosmology. The solid an- 
gle of the lightcone is 75 x 15 degrees which is split into three 
75x5 degree slices. The 2QZ Survey consists of 2 such slices. 
Ideally we would like many more than 3 slices but due to 
the large volume of the survey (4.2xl0 9 /i~ 3 Mpc 3 for z < 2.2 
and f2 m (0)=0.3, f2A(0)=0.7) this is currently not possible. 



2.3 Construction of the Mock Catalogues 

To create realistic mock catalogues, there are a number of 
steps which must be carried out. First the mass particles 
must be biased in order to match the predicted clustering 
amplitude of 2QZ QSOs. Initial analysis of the correlation 
function, assuming a cosmology with f2 m (0) = 0.3, S7a(0) = 
0.7, found that the correlation function could be well ap- 
proximated by a power law of the form £(s) = (s/so)~ 7 , 
with 7 rs 1.7 and sp ~ 6.0/t~ 1 Mpc, at all epoc hs, consistent 

Shanks (1996). See proom et al 



with the work of 



Croom 



(2001a) for the correlation function of the 10k catalogue of 



2QZ QSOs. 

The biasing prescriptio n used in th i s pap er is similar 
to the method described in Hoyle et al. (199SJ) and follows 
method 2 as described by Cole et al, (p.998[). The bias prob- 



ability is based on the density field at the epoch at which 



* The Hubble Volume simulations were performed by the "Virgo 
consortium for cosmological simulations". This is an interna- 
tional collaboration involving universities in the UK, Germany 
and Canada. The members of this consortium are: J. Colberg, H. 
Couchman, G. Efstathiou, C. Frenk (PI), A. Jenkins, A. Nelson, 
J. Peacock, F. Pearce, P. Thomas, and S. White. G. Evrard is 
an associate member. The Hubble Volume simulation was car- 
ried out on the Cray-T3E at the Max-Planck Rechen Zentrum in 
Garching. 



particles are selected rather than from the initial density 
field. The bias is approximately scale independent on the 
scales we probe here. A full d escrip tion of the biasin g proc e- 
dure can be found in Hoyle (2000J) and Hoyle et al. (2001). 

We also have to match the radial selection function of 
the QSOs. We fit a polynomial to the N(z) distribution of 
the QSOs with redshifts in the range 0.3 < z < 2.2 to ob- 
tain a probability distribution function and randomly select 
particles using this until there are «12,500 mock QSOs on 
each slice. 

To mimic the two independent slices of the 2QZ, we 
measure £(cr, n) from the two outer slices and average the 
results together. 



3 MEASURING £(cr,7r) 

There are many different ways in which the clustering per- 
pendicular (a) and parallel (7r) to the line of sight can be 
defined. Here, we define tv = \s2 — Si| and a — (si + S2)9/2, 
where si and S2N are the distances to two QSOs and 9 is 
the angle between them, measured from the position of an 
observer. 

£(a, n) is then estimated in much the same way as the 
two-point correlation function. A catalogue of unclustered 
points, that have the same radial selection function and an- 
gular mask as the data but which contains many more points 
than the data catalogue, is used to estimate the effective 
volume of each bin. The DD(a,n), DR(a,n) and RR(a,n) 
(where D stands for data or mock QSO and R stands for ran- 
dom) coun ts in each -k and a bin are found and the Hamilton 
estimator ( Hamilton! 1993 ) 

DD{a,n)RR(a,n) 



£(cr,7r) = 



DR(o-,tt) 



(1) 



is used to find £(<r, n). Due to the sparsity of QSOs 
in the 2QZ Survey, we use bins of <51og(7r//i -1 Mpc) = 
51og((j//i _1 Mpc) = 0.2. The sparsity of QSOs also means 
that Poisson errors are a re asonable esti mate of the error on 
scales less than 40/i _1 Mpc ( jHoylcj |200c| ) . 



A£(cr,7r) = (l + £(cr,7r)) 



DD(a, 7r)' 



(2) 



The factor of two is needed as we require the number of 
independent pairs in each bin. 



4 MODELING £(<r, n) 

4.1 Definitions 

There are many terms used in this analysis that can be easily 
confused. We define their meaning here and stick to these 
conventions throughout. 

• Underlying cosmology - this is the (unknown) cosmol- 
ogy of the Universe 



T We adopt the convention that s refers to apparent distances in 
redshift space and r refers to distances in real space through out 
this paper. 
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• The simulation cosmology - the known cosmology 
(fi m (0)=0.3, n A (0)=0.7) of the Hubble Volume Simulation 

• Assumed cosmology - the cosmology used when mea- 
suring the two-point correlation function and f(cr, n) from 
the 2QZ Survey or the Hubble Volume simulation. Models of 
£(<7, 7r), which are discussed later, also have to be calculated 
in the same assumed cosmology as the data. We consider two 
possibilities for the assumed cosmology, fi m (0) = l, fi A (0)=0 
(EdS) or Q m (0)=0.3, fi A (0)=0.7 (the A cosmology), when 
fitting the models to the simulation to show that the results 
are not sensitive to this choice. 

• Test cosmology - the cosmology used to generate the 
model predictions for £(<r, 7r) which are then translated into 
the assumed cosmology. When the test cosmology matches 
the underlying (or simulation) cosmology, the distortions in- 
troduced into the clustering pattern should be the same in 
the model and in the data. The model should then provide a 
good fit to the data, providing the redshift space distortions 
have been properly accounted for. If the test cosmology is 
incorrect then the model should not fit the data. This po- 
tentially allows the underlying cosmology to be determined. 



4.2 Requirements of the Model 

One of the aims of this paper is to test whether the cosmo- 
logical parameters, f2 m (0), S1 A (0) and also (3qso{z) will be 
constrained from the 2QZ Survey. /3 g (0) has been measured 
from galaxy redshift surveys by comparing the ze roth and 
second order moments of the correlation function ( Ratcliffc 



et al 



1998c). However, the effects of cosmology add an extra 



distortion to the clustering pattern of objects with a wide 
range of redshifts. This alters the moments of the clustering 
too, such that the effects of redshift space and cosmological 
distortions are difficult to disentangle. 

Instead, we test if cosmology and Pqso (z) will be con- 
strained by comparing £(<r, n) measured from the mock cat- 
alogues to models of £(ct, n). The actual procedure is de- 
scribed in more detail in section H. The idea is that when 
the values of fi m (0), fi A (0) and j3qso{z) used to calculate 
the models are the same as the underlying values, the model 
will match the mock catalogues, allowing the cosmology to 
be determined. 

There are several assumptions that go into the models. 
One of the assumptions is that we are comparing the models 
to the data on linear scales so the effects of non-linear clus- 
tering can be ignored in the models. A second assumption 
we make is that the bias factor, b, is independent of scale. 
This is the case, by design, in the mock catalogues but this 
may only be an approximation on small scales in the 2QZ 
Survey. 

The lit erature contains many e x amples of models for 



Ratcliffe et al 



Matsubara fc Suto 1993, Popowski et al 



1998 



1998c[) and the power spectrum, measured 



parallel and perpendicular to the lin e of sight, P{k\\,k± r , 
(Ballinger, Peacock & Heavens 1996). We consider various 



aspects of these models as the models must account for the 
effects of cosmology and the redshift space distortions. We 
describe the method for including the cosmology below and 
describe a model for the effects of redshift space distortions 
in section 4.4. 



4.3 Effect of A on Clustering Anisotropy 

If the assumed cosmology is different from the underlying 
cosmology of the Universe (or the simulation), distortions, 
different from those caused by peculiar velocities, are intro- 
duced into the clustering pattern. This is because the radial 
and perpendicu lar directions are affected by c o smolo gy in 
different ways. Ballinger, Peacock & Heavens (1996) 



out- 



line how the cosmology scales the power spectrum split into 
its parallel and perpendicular components. Cosmology sim- 
ilarly affects the correlation function in the linear regime 
according to: 



U 



(<7,7r) =£ ta 



(3) 



i.e. the correlation function in the assumed cosmology is 
the same as the correlation function in the test cosmology 
provided a and it are scaled according to 

' a aB » ( 4 ) 



and 



f± 

TV 
J 



B t 

7T A a 



(•») 



where the subscripts a and t refer to the assumed cosmology 
and the test cosmology respectively. For a flat Universe, A 
and B axe defined as 

A =4- (6) 



H o v /fi A (0) + O m (0)(l + , 



and 



dz 



(7) 



B= _£_ 

H °Jo \A2 A (0)-rfi m (0)(l + ,-'V ; 

These can be calculated for open universes too -- see 
Ballinger, Peacock & Heavens (1996) for the full definition 



— however we only consider flat universes. Both A and B 
are calculated at the median redshift of the survey, which in 
this case is z — 1.4. By comparing the models to the sim- 
ulation, we find that this is an adequate approximation to 
make (see Figure III). 

4.4 Redshift Space Distortion Model 



Following |Peebles| ( ft980| ) and ftatcliffe et al.| ( |1998c[ ) we de- 
fine the correlation function, parallel and perpendicular to 
the line of sight by 



l+£(a,7r)= [1 + ar)]g(r,w)dw :i 



(8) 



where f (r) is the real space QSO or mock QSO correlation 
function, free from the effects of redshift space distortions. 
w = vi — Vj where v is the peculiar velocity of a QSO after 
subtracting off the Hubble flow and r = ri — rj. g(r, w) is 
the distribution function of w for QSO's separated by r. 
Here r 2 = o 2 + r 2 , where r z — ty — w z /H and w z is the 
component of w parallel to the line of sight, denoted by z 
for simplicity. 

If it can be assumed that g(r,w) is a slowly varying 
function of r then g(r, w) = <?(w). Then equation w can be 
simplified and becomes 



l + f(a,7r) 



[1 + Z(r)]f (w z )dw z 



(9) 
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where 



f(w z ) = / Aw x I dw y g(w) 



(10) 



A simple streaming model for the bulk motions of galax- 
ies can be included by writing g(r, w) = g[w — rv(r)] where 
v(r) is the mean relative radial velocity of galaxies separated 
by r. If the approximation for g(r, w) in equation hoi is made 
again, equation H becomes 



i+e(*,7r) 



[l + $(r)]f[w s -v(r z )]dwz 



(11) 



Models for the effects of small scale peculiar velocities 
and the bulk motions of galaxies or QSOs are required if 
£(<7, 7r) is to be accurately described. 



4-4-1 Fingers of God 

The distribution function of the small scale peculiar veloc- 
ities has been previou sly modeled as an exp onential or a 
Gaussian distribution. Ratcliffe et al (1998c) found that a 
substantially better fit to peculiar velocities in TV-body sim- 
ulations was found if an exponential model was used and 
that is what we adopt here: 



f(u>z 



1 



V2<W 2 Z >V2 



exp 



-y/2 



< wl > 1 / 2 



(12) 



where < w 2 > x / 2 is the rms line of sight pairwise velocity 
dispersion. 

Errors in the measurements of the QSO redshifts may 
also affect the measurement of £(a, n) on small scales. These 
errors may lead to the peculiar velocities being better mod- 
eled by a G aussian distribution. We discuss this further in 
Section 8.1. 



4-4.2 Bulk Flows 

The model that we use for the bulk motions depends on the 
cosmology of the Un iverse, the cl u sterin g of QSOs and the 
QSO bias. Following |Hale-Sutton] ( |L99C| ) we set 



«(>"*) = -; 



'7 



-CL 



l (0)°-°H(z)r M 



m 



b 2 + £(r) 



(13) 



r z is defined above, /3 and fi m (0) are free pa rameters an d 
the bias, b, is calculated using b — £7 m (0)°' 6 //3 (Kaisei 1987). 



Clearly in the case of the QSOs, b is the QSO bias and b, Cl m 
and j3 are defined at the average QSO redshift. The value 
of H a in the simulation is 70 km s _1 Mpc -1 , which we will 
assume when we analyse the data too, and the evolution of 
Hubble's constant with redshift and fiat cosmology is given 
by 

H(z) = H {n m (0)(l + zf + fi A (0)}°' 5 . (14) 

We self-consistently determine the real space correlation 
function from the redshift space correlation function. This is 
discussed in section 0. The real space correlation function is 
not approximated by a power law but a value of 7 is still re- 
quired in equation |l3|. A value of 7=1.7 best fits the redshift 
s pace correla tion function over the range 4 < s < 40/i~ 1 Mpc 
( Hoyle 2000 ) so this is the value used in equation [13. 



An alternative method for modeling the e f fects o f in 
fall on £(cr, 7r) is outlined in Matsubara fc Sutc ( 1996). Es - 
sentially, they generalise the formula of Hamilton (1995 



to define a relation between the redshift space correlation 
function at z / and the real space correlation function. 
This model includes the effects of infall but not small scale 
peculiar velocities. These could be included by convolving 
the redshift space correlation function with an exponential 
model for small scale peculiar velocities along the 7r direc- 
tion. 



4.5 Comparison with the Hubble Volume 
Simulation 

Does the model reproduce £(c, 71") measured from the mock 
QSO catalogues? The left hand plots in Figure Fy show f (<r, 7r) 
measured from the mock catalogues. We show the fractional 
errors on £(er, 7r) in the middle plots — the lighter the shad- 
ing, the smaller the errors. The areas of the £(c, 7r) diagram 
where the fractional errors are the smallest are the areas 
where most of the differentiation between different models 
of £(<r, 7r) can be made. The right hand plots show the model 
£(<j, 7r). Three different cosmologies have been assumed for 
the simulation and the models, f2 m (0) = l, $1a(0)=0 (shown 
in panels a, d and g), fi m (0)=0.3, f2A(0)=0.7 (panels b, e and 
h), which is also the simulation cosmology, and f2 m (0) = 0, 
Oa(0) = 1 (panels c, f and i). The bold, solid lines show 
£=0.1, the dot dashed bold lines show £=1. The solid lines 
increase from £=0.1 in steps of 0.1, the dashed lines decrease 
from £ = 0.1 in steps of 0.01 and the dot-dashed lines in- 
crease from £=1 in steps of 1. The values of the velocity 
dispersion (< w\ > 1/2 =400 km s" 1 ), /3qso(z) (=0.36) and 
test cosmology in the models are the values measured from 
the Hubble Volume simulation at the average redshift. 

The measurement of £(<r, tv) from the Hubble Volume 
simulation is quite noisy. However, particularly if the bold, 
solid line (£ = 0.1) is considered, the effects of the different 
assumed cosmologies can be seen. This line in panel a) is 
far more squashed in the -k direction than in panel b) and 
elongation is seen in panel c) as compared to panel b). This 
is seen in the models too. The solid line in panel d) appears 
more squashed in the 7r direction than the lines in panel e) 
and the lines in panel f) appear more elongated in the 7r di- 
rection than the lines in panel e). For each choice of assumed 
cosmology, the models match the simulation to within ~ la. 



5 DETERMINING THE UNDERLYING 
(SIMULATION) COSMOLOGY 

We wish to use geometric distortions in the measurement 
of £(<t, 7r) to constrain the cosmological parameters fi m and 
Qa- We apply the method to QSOs as it is only at high 
redshifts that the geometric distortions become apparent. 
Current, wide-angled galaxy redshift surveys only probe the 
universe out to z ~ 0.1 — 0.2 where an incorrect choice of 
cosmology leads to a 10% uncertainty in separations for flat 
cosmologies. This produces only a small distortion in the 
clustering pattern which is very hard to seperate from the 
effects of bulk flows. However, if clustering is measured over 
a wider range of redshifts, the effects of assuming the wrong 
cosmology become large, as can be seen in Figure hi 
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Figure 1. The left hand plots show £(tr, 7r) measured from the mock catalogues. The assumed cosmology is C m (0)=l, f2A(0)=0 (top 
row), £7 m (0)=0.3, Ca(0)=0.7 (middle row) and O m (0) = 0, Oa(0) = 1 (bottom row). Each strip contains 12,500 mock QSOs and the 
average over the two outer strips is shown. The solid, bold contour shows £ = 0.1, with the other solid lines increasing logarithmically 
from it with 0.1 interval and the bold, dot-dashed contour shows £ = 1. See the text for the value of the other lines. The centre plots 
show the fractional errors on §(cr, it). The white area shows where £(<r, 7r) is measured to better than 40% and with each contour the 
fractional error increase by 20%. The right hand plots show models of £(<r, it). The values of the velocity dispersion, /3qso(2) and the 
test cosmology are taken from the Hubble Volume simulation. 



We constrain cosmology by comparing £(cr, n) measured 
from the data, or here the Hubble Volume, to our models. 
There are a number of free parameters in the models, the test 
cosmology, Pqso(z) and the small scale peculiar velocities, 
< w 2 > 1 ' 2 . We fix the assumed cosmology to be either EdS 
or the A cosmology to show that the results are independent 
of this choice. 

As we compare the models and the simulation on scales 
greater than 4/i _1 Mpc in both the 7r and a direction, we 
find that the modeling of £(<r, n) is fairly insensitive to the 
value of the small scale velocity dispersion. This is shown 
in Figure 0. Here we show two different sets of models, 
the top plot shows models with test cosmology f2 m (0) = l, 



Aa(0)=0, /3qso(z)=0.36 and the A cosmology assumed. 
The solid lines have < w\ > ' 2 — 50 km s~ x and the 
dashed lines have < w\ y 1 ' 2 — 1000 km s _1 . The lower plot 
shows models with test cosmology £7 m (0)=0.3, fiA(0)=0.7, 
/3qso(z)=0.36 and the A cosmology assumed. Again the 
solid lines have < w'j. y 1 ' 2 — 50 km s _1 and the dashed lines 
have < wl > 1 ' 2 = 1000 km s _1 . Considering the wide range 
of velocities covered by these two models, there is little dif- 
ference between them, apart from on the smallest a scales 
where the errors on £(<r, n) from the mock catalogues are 
large anyway (see Figure hi). The models with both choices 
of < w'j. y 1 ' 2 still match the simulation £(cr, w) to within la. 
This is because most of the differentiation between different 



Constraining Cosmology from Redshift Space Distortions via £(a, it) 7 



models of £(a, it) occurs on larger scales. We therefore fix 
< w\ y 1 ' 2 — 400 km s _1 in the models when we are finding 
the Hubble Volume cosmology. 

As the effects of small scale peculiar velocities have little 
impact on models of £(a, 71") on scales larger than 4ft _1 Mpc, 
they should also have little effect on the correlation func- 
tion £(s). We have tested this on the Hubble Volume simu- 
lation by comparing the real and redshift space correlation 
functions of the mass. We find that on scales greater than 
~ 3ft Mpc the small scale peculiar velocities have negligi- 
ble effect on the shape of the redshift space mass correla- 
tion function. It is therefore possible to find the real space 
correlation function self-consistently from a redshift space 
co rrelation func tion for each test value of j3 via the formula 
of |Kaiser| (|l987|) 



«r) 



«*) 



{^ + ¥ + ¥ 2 Y 



(15) 



where r and s denote the correlation function in real and 
redshift space respectively and /3 = fi m (0) 0,6 /&. This is mea- 
sured in the assumed cosmology but, to create the model 
£(<7, 7r), we need the real space correlation function mea- 
sured in the test cosmology. First we translate the redshift 
space correlation function, £(s), in the test cosmology to the 
assumed cosmology, £(s') by scaling s to s' using 

2 



2 . 2 

a +r z 



and 



W/x) 2 + (^//||) 2 



(16) 



(17) 



then we apply equation h_5l The terms f± and f\\ are defined 
equatio n r& and H as well as in Ballinger, Peacock & Heav- 



ens (19 )6). The accuracy of converting the two-point corre- 
lation function from one cosmology into another is demon- 
strated in Figure H. The symbols show the redshift space 
correlation function measured from the Hubble Volume with 
EdS assumed (triangles) and the A cosmology assumed (cir- 
cles). The line shows the correlation function measured in 
the EdS cosmology but translated into the A cosmology. 
There is excellent agreement between the two over a wide 
range of scales. This is a further consistency check that mea- 
suring fx and /|| at the average redshift of the Survey is an 
adequate approximation to make. An alternative but more 
time consuming method would be to measure the redshift 
space correlation function in each of the test cosmologies di- 
rectly and then apply equation [la to obtain the real space 
correlation function in the test cosmology. 

We now just have two free parameters, the test cosmol- 
ogy and the value of /3qso{z)- We run a grid of models and 
find which choice of test cosmology and /3qso(z) give the 
best fit to the data. 



5.1 The Fitting Procedure 

The fitting procedure that we adopt to find which test cos- 
mology matches the simulation cosmology is as follows: 

1) Pick an assumed cosmology (here either EdS or the A 
cosmology) 

2) Calculate £(s) and £(<r, 7r) from the data (the 2QZ Survey 
or the mock catalogues) using the assumed cosmology. 

3) Pick a value for the test /3qso(z) at the average redshift 
of the survey. 




Figure 2. Two sets of models with different values of the small 
scale peculiar velocities. The top plot shows models with test cos- 
mology S7 m (0)=l., Ca(0)=0., [3qso(z)=0.36 and an assumed A 
cosmology. The solid lines have < w 2 > 1 / 2 = 50 km s _1 and 
the dashed lines have < w 2 > 1 / 2 = 1000 km s — 1 . The lower 
plot shows models with test cosmology n m (0)=0.3, nA(0)=0.7, 
/3q5q(z)=0.36 and an assumed A cosmology. Again the solid 
lines have < w 2 > 1 / 2 = 50 km s — 1 and the dashed lines have 
< w 2 > 1 / 2 = 1000 km s" 1 . Apart from on the smallest scales, 
there is little difference between the two models considering the 
large difference between the two velocities and the error contours 
shown in Figure hi 



4) Pick a value of the present day test f2 m (0) and calculate 
n m (z) at 2=1.4 via fi m (z) = fi m (0)(l + z) 3 /[fim(0)(l + z) 3 + 

n A (o)(o)]. 

5) Calculate the bias b QSO {z) = 0. m iQf- e {z) / P QS o{z). 

6) Translate the redshift space correlation function, £(s), 
measured in the assumed cosmology to a real space correla- 
tion function, £(»"'), in the test cosmology using equation |l7| 
and|I|. 

7) Generate the model £(<r, 7r). 

8) Translate the model ^(a, n) from the test cosmology into 
the assumed cosmology using equations H and m 

9) Calculate how well the model £(<r, n) fits the data £(cr, n) 
via the x 2 statistic, using the Poisson errors from the data 
£(a, 7r) measured in the assumed cosmology. 

10) Go back to 3) using a different test cosmology and 
Pqso{z). 

When the parameters /3qso(z) and fi m (0) nratch those 
of the underlying cosmology (or in this case the simulation 
cosmology), the value of \ 2 should be minimised. We fit the 
model £(<t, -k) to the Hubble Volume ^(a, -k) over the range 
of scales 4< a, ty <40ft _1 Mpc. This is to ensure that any 
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Figure 3. The redshift space correlation function of the mock 
catalogues estimated assuming the A cosmology (circles) and EdS 
(triangles). The line shows the correlation function measured as- 
suming EdS corrected to the A cosmology with the relation be- 
tween s and s' given in equation rL 71. This line agrees well with the 
circles over the range of scale 3< r <50h~ 1 Mpc, showing that a 
correlation function measured in one cosmology can be translated 
into another. 



non-linear effects are small and that the errors on £(cr, n) do 
not dominate the actual value of £(<x, 7r). 

The error contours, shown in Figure W are found as fol- 
lows. The number of degrees of freedom, v, is the number of 
bins in which the model £(cr, n) is fitted to the data £(a, n) 
minus the number of free parameters (this is 2, f2 m (0) and 
Pqso{z))- We have checked the error contours using Monte 
Carlo simulations and find the two methods to give compa- 
rable results. 



6 PREDICTED RESULTS FOR THE 2QZ 
SURVEY 

The results of fitting the £(a, tt) models to £(a, n) measured 
from the mock catalogues are given in Figure H. The up- 
per panel (a) shows the results when the assumed cosmol- 
ogy is EdS and the lower panel (b) is for the A cosmol- 
ogy. We show the results for two different assumed cosmolo- 
gies to show they are not dependent on this choice. In both 
cases, the simulation cosmology, £7 m (0)=0.3, fiA(0)=0.7 and 
Pqso(z) — 0.36 is contained within the la contour. How- 
ever, the results do not place a strong constraint on fi m (0) 
as there is a degeneracy between Q m (0) and Pqso{z). How- 
ever, the value of Pqso{z) with z — 1.4 should be fairly well 
constrained from the 2QZ Survey. This does at least place a 
joint constraint on cosmology and the QSO-mass bias. 

The reason for this degeneracy is that the errors 
on £(a, tv) from the mock catalogues are fairly large. 
Shown in Figure U are two different $,(a, n) models. One 
has On(0)=0.3, Ov(0)=0.7 and a value of f3 QSO (z)=0.35 
(dashed lines) and another model has fi m (0) = l, Qa(0)=0 
and [3qso(z)=0.5 (solid lines). All the other parameters (the 
assumed A cosmology and the small scale peculiar velocities 
of < w'i > 1 ' 2 =400 km s _1 etc) are the same in each model. 





Figure 4. Fitting £(<r, 7r) from the Hubble Volume. The top panel 
(a) has EdS as the assumed cosmology and the lower panel (b) has 
the A cosmology assumed. In both cases the simulation cosmology, 
fim(0)=0.3, Oa(0)=0.7 and /3qsq(z)=0.36 is contained within 
the la contour, although little constraint on f2 m (0) is possible 
from the mock 2QZ Survey catalogues. 



These models both provide a good fit to the mock catalogues 
from the Hubble Volume (see Figure vt) which is not surpris- 
ing as the two models have similar shape and amplitude. 
Very small errors on £(a, n) are needed to detect the small 
differences between these models. 

The degeneracy between Si m (0) and /3qso(z) seen in 
Figure does not lie in the direction that is perhaps in- 
tuitively expected. Consider the two assumed cosmologies 
adopted here, the two test cases shown in Figure q and the 
simulation cosmology fi m (0)=0.3, £Ia(0)=0.7. If EdS is as- 
sumed then £(cr, 7r) from the simulation is squashed (seen 
in Figure Dla)), as is the model £(<r, tv) with test cosmology 
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Figure 5. The dashed lines show a model with test cosmol- 
ogy Q m (0)=0.3 and C A (0)=0.7 with a value of /3 QSO (z)=0.35. 
The solid lines show a model with test cosmology f2 m (0)=l and 
/3qso(S)=0.5. The difference between the two models is small. 
The A cosmology is assumed in both cases. 



fi m (0)=0.3, Oa(0)=0.7, due to the simulation and test cos- 
mologies being different to the assumed cosmology. £(cr, n) 
using fi m (0) = l, f2A(0)=0 as the test cosmology is not com- 
pressed as the test cosmology matches the assumed cos- 
mology. A higher value of Pqso (z) is therefore required to 
squash the model £(a, n) to match the simulation. 

If the A cosmology is assumed, then no squashing due 
to cosmology occurs on the simulation £ (a, tv) (again seen 
in Figure IlKb)) or on the model £(<t, 7r) with f2 m (0)=0.3, 
Qa(0)=0.7 as the test cosmology as in this case the simu- 
lation and test cosmologies match the assumed cosmology. 
However, if fi m (0) = l, Oa(0)=0 is the test cosmology, then 
the model £(a, tt) is elongated due to the different test and 
assumed cosmologies. Therefore, a larger value of /3qso{z) 
is again required to compensate for the elongation, allowing 
the model to match the simulation £(cr, 7r). This is partly 
shown in Figure H as here we have models with the A 
cosmology assumed and /3qso(z)=0-35 (although the small 
scale peculiar velocities are 50 and 1000 km s _1 rather than 
400 km s _1 ). The top plot has test cosmology fi m (0) = l, 
?2a(0)=0 and the lower plot has test cosmology f2 m (0)=0.3, 
^a(0)=0.7. Elongation is clearly seen in the top plot as 
compared to the bottom plot, showing that a higher value 
of j3qso{z) is required for the model with test cosmology 
n m (0) = l, Oa(0)=0 if the model is to match the simulation. 

The results for determining the cosmology from the 
mock catalogues suggest that the cosmological parameters, 
fi m (0) and Oa(0), will not be well constrained from the 2QZ 
Survey. However, the fits are only made to the mock 2QZ 
catalogues and there may be ways in which the constraints 
from the final 2QZ may turn out better than these results 
suggest. 

The corre l ation function from the 2QZ 10K Catalogue 
(Croom et al. 2001a) has a powerlaw form over a slightly 



wider range of scales than the correlation function from the 
Hubble Volume simulations, which turns over to lower am- 
plitude at small and large scales. Therefore, on small and 
large scales, the correlation function from 2QZ has a slightly 
higher amplitude than that from the Hubble Volume. This is 
mainly because the biasing scheme was chosen to match the 
clustering of the 2QZ correlation function on scales around 
10/i _1 Mpc but we also preserved the shape of the dark mat- 
ter correlation function. The higher amplitude in the 2QZ 
correlation function on small and large scales helps in two 
ways. The extra large scale power may allow us to fit the 
models out to larger scales, increasing the number of bins 
over which we can compare the data and models. Also, the 
number of pairs found in the 2QZ Survey on small scales 
will be larger than the number found in the mock cata- 
logues. This reduces the errors on small scales, hopefully 
giving tighter constraints on the models. 



7 FURTHER CONSTRAINTS ON 
COSMOLOGY AND /3qso(z) 

As shown in section n, the cosmological parameters, fi m (0) 
and I3qso{z), will probably not be uniquely measured from 
the 2QZ Survey using this technique. However, there are 
other ways in which the cosmological parameters can be es- 
timated from the 2QZ Survey and results can be combined 
to obtain stronger constraints on cosmology. The basic idea 
developed in this section is that the f2 m : /3qso(z) degen- 
erate set obtained from QSO clustering evolution is differ- 
ent from the f2 m : I3qso(z) degenerate set obtained from 
analysing QSO space distortions and that by using these 
two constraints in combination, the degeneracy in these two 
parameters may be lifted. 

The amplitude of QSO clustering can be measured at 
z — 1.4 assuming any cosmology. In order to calculate the 
value of /3qso(1-4), we need to know the clustering ampli- 
tude of the mass at this redshift. To obtain this, we need to 
know the amplitude of the clustering of the mass at z — 0, 
then linear perturbation theory allows us to trace the evo- 
lution of the mass clustering. 

However, the clustering of the mass at z — is not 
known. Here we use information from galaxy redshift surveys 
in order to determine the mass correlation function at z = 
0. The mass correlation function can be determined if the 
galaxy correlation function and /3 g are known, assuming that 
bias is independent of scale. 

Many measureme nts of /3 R at z = have appea r ed in the 
literature, for e xample Tadros fe Efstathiovj ( 1996 ^_ Ratcliffe 



et al 



1998c|) [Hoyle ct al.| ( |1999| ) and putram, Hoyle and 



Shanks| (|200C). The latest measurement of /3 S co mes from 



the 2dF Galaxy Redshift Survey and is 0.43T0.07 (Peacock 



et al 



2001). For each cosmology (again we only consider flat 



cosmologies) the value of the galaxy-mass bias can be found 
from /3 g , which in turn gives the value of the mass correlation 
function if the galaxy correlation function at z = is known, 
as shown below. 

Rather than determining the value of the two-point 
correlation function at one particular point, we use the 
less noisy volume averaged two-point correlation func- 
tion, £|. We assume that the galaxy correlation function 
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can be approximated by a power law of the form £| 
(s/6h~ 1 Mpc)~ 17 . |g is then found via 



£ 3 = 



^U^is/ey^ds 



f 

Jo 



47rs 2 ds 



(18) 



By integrating out to 20/i _1 Mpc, the non-linear effects on 
the volume averaged correlation function should be small. 
The power law approximation is in very good agreement 
with early results from the 2dF Galaxy Redshift Survey 
two-point correlation function out to 20/i~ 1 Mpc (P. Nor- 
berg, private communication). Once the 2dF Galaxy Red- 
shift Survey is finished, there will be no need to make this 
approximation. 

£qso i s found from the Hubble Volume by finding the 
number of DD, DR and RR pairs that have separations in 
the range 0-20/i -1 Mpc and forming the two-point correla- 
tion using the estimator of Hamilton. These values agree well 
with the values obtained from the 2QZ 10K Catalogue as 
the mock correlation function was biased to approximately 
match early results from the data at ~10/i _1 Mpc. 

For each cosmology, the bias between the galaxies and 
the mass at z=0, b gp (0), can be found from 



Mo) 



»°; 6 (Q) 
&(o) ' 



(19) 



The real space galaxy correlation function can be deter- 
mined from the redshift space galaxy correlation function 



£ r = 






[l + |/3 g (0) + l/3|(0)] : 



(20) 



where the superscripts r and s indicate real and redshift 
space respectively. The real space mass correlation function 
at z=Q can now be found as 



£(0) 



ML 

MO)' 



(21) 



The real space mass correlation function evolves according 
to linear theory such that 



£(*) 



g(Q) 
G( z y 



(22) 



where G(z) is the growth factor, which depends on cos- 
molog y, found from the formula of Carroll, Press & Turner 
(h.992h. When fi m (0)=l, G{z) = (1 + z). 



We now relate the correlation function of the mass at 
z measured in real space to the amplitude of the QSO clus- 
tering at z, measured in redshift space as we wish to know 
Pqso(z) as a function of f2 m (0). First we calculate fi m (,z) 
using 



fim(z) = 



n m (o)(i + z ) 



fi m (0)(l + z) 3 + O\(0)' 



(23) 



This is valid for flat cosmologies only. (3qso( z ) is then given 
by 



Pqso(z) - 



n m («) 



but as yet b qp (z) is unknown. b qp (z) is defined as 
,2 , n £qso( z ) 

Kpiz) = "W 



(24) 



(25) 



substituting the above equations into equation h4 gives 

<iQSO\ z ) 

£qso( z ) can be found via 

* , n _ £qso( z ) 

<iQSO( z 



(i + y Q so{z) + ^ QSO (z)y 

Substituting this into equation |2(j, gives 
Pqso( z ) =H m (z) 



(26) 



(27) 



1 ' 2 Tr^TT[ 1 +|^oW+^QSo(^)],(28) 



*QSO\ 

a quadratic in /3q SO { z ) which can easily be solved, allowing 
Pqso{ z ) for any il m (0) to be found. We substitute z = 1.4 
to find the value of /3qso(z) at the average redshift of the 
survey. 

The errors on /3qso( z ) are found in the standard way, 
i.e. by differentiating (3qso( z ) with respect to all the vari- 
ables that contribute to the error and summing the compo- 
nents of the error in quadrature. In this case there are errors 
on /3 9 (0), £qsq{ z ) and g(0). The error on (3 g (0) is ± 0.07, 
the error on £qso( z ) is the Poisson error found from the 
total DD counts. The error on £ s is the least well known as 
the power law approximation is used. We assume an error 
of around 20%, which is larger than the error expected from 
the final 2dF Galaxy Redshift Survey. 

7.1 Model Results 

A comparison between the constraint on f2 m (0) and Pqso( z ) 
found here and the constraint found from fitting models 
to the full 2QZ Hubble Volume mock catalogue £,(a, n) is 
shown in Figure n. The contours are shown as in Figure 
H and the solid lines show the la and 2a constraints from 
the above method. We show the comparison for the two as- 
sumed cosmologies. It can be seen that f2 m (0)=0, fiA(0) = l 
and £7 m (0) = l, J1a(0)=0 are ruled out with ~ 2a confidence 
when the results are combined. The value of f2 m (0) is now 
constrained to be ~0.3l ' 15 , by considering where the la 
contours from the two estimates overlap. 

The constraints on /3qso( z ) are even stronger than 
the constraints on cosmology. /3qso( z ) can be constrained 
from the £(ct, 7r) contours alone but by combining the er- 
rors, /3qso( z ) can be measured to an accuracy of ~ ±0.1 
which places joint constraints on cosmology and the QSO- 
mass bias at the very least. This suggests that when 2QZ 
is finished, a powerful method for constraining cosmological 
parameters will be available. 



8 THE 10K CATALOGUE 

Our overall aim here is to apply th e analysis described above 
to the 10k catalogue of 2QZ QSOs flCroom et alJbOOlbh. The 



catalogue contains 10681 QSOs and is publically available 
at www.2dfquasar.org. This catalogue contains the most 
spectroscopically complete fields (i.e. fields in which more 
than 85 per cent of objects have been identified) that were 
observed prior to November 2000. 

As the survey is not finished, it currently has a patchy 
angular selection function. A 1' by 1' completeness map is 
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Figure 6. Combining predicted results from measuring £(<r, it) 
from the Hubble Volume mock catalogues of the full 2QZ Survey 
(contour plot as in Figure W) with results from the evolution of 
clustering. Plot (a) shows the contours with EdS assumed and plot 
(b) shows the contours with the A cosmology assumed. The lines 
show the 1 and 2<r contours predicted from the method assuming 
/3 g al- The dashed lines show the la joint contours. 



constructed for both the north and south galactic cap (NGC 
and SGC hereafter) regions to allow us to construct a ran- 
dom catalogue with the same angular selection function. We 
also take into account the effects of galactic extinction in the 
co nstruction of t he ran dom catalogue. Fu l l deta ils are given 
in |Croom ct ai] ( fc001a| ) and JHoyle et alj (feoOlj ), 

There is an uncertainty in the measurement of the QSO 
redshifts due to the resolution and signal to noise of the 
2dF spectra. Redshifts can currently only be measured t o 
an accuracy of at best a{z) = 0.0035 (Croom et al. 2001b). 
This means that as well as the intrinsic small scale 'Finger of 




Figure 7. Figure (a) shows the constraints possible when we 
apply the 10k masks to the Hubble Volume cosmology. We assume 
the A cosmology. The patchy angular selection function does not 
seem to bias the results. In (b) we apply the same analysis to 
the 10k masked Hubble Volume but we have added a 5h~ 1 Mpc 
Gaussian dispersion to the positions of the QSOs to mimic the 
redshift errors present in the 10k 2QZ sample. The dashed lines 
show the joint ltr contours. In both cases we are still able to 
recover the cosmology of the Hubble Volume. 



God' velocity dispersion, there is also an apparent velocity 
dispersion of around 600 km s _1 . If QSOs have a similar 
velocity dispersion as galaxies of ~ 500 km s _1 (Ratcliffe 
et al. 1998c) then the observed velocity dispersion would be 
around 700-800 km s" 1 . 
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Figure 8. The solid lines show £(<r, it) measured from the 10k data catalogue assuming the EdS cosmology (a) and the A cosmology 
(b). The dashed lines show models. The models have the best values of f2 m (0) and /3qso(^) found by combining the two techniques of 
constraining cosmology. When we assume EdS the values are f2 m (0)=0.20, /3qso=0.35, when we assume the A cosmology the values are 



Q m (0)=0.3, /3qso=0.35. We set < w\ > 1 /2 = goo km s" 1 in the models. 



8.1 10k Hubble Volume Mock Catalogue Analysis 

To test the effect of the current window function on the 
fitting of £(a, it), we have again used the Hubble Volume. 
We imprint each of the completeness maps (NGC and SGC) 
onto one of the Hubble Volume slices, measure £(a, n) and 
combine the results as we do with the data. Due to the 
smaller number of QSOs in the current catalogue, the mea- 
surement of £,(a, it) is more noisy than that from the fi- 
nal catalogue will be. Therefore, rather than binning in 
<5rog(er, it) — 0.2, as was done for the analysis of the 25k mock 
catalogues, we use slightly larger bins of <51og(cr, it) — 0.3 for 
the 10k mock and 2QZ analysis. This increases the signal 
within each bin but reduces the number of bins available for 
comparison with the models. As the uncertainty in the red- 
shifts increases the small scale peculiar velocities, we restrict 
our fitting range to 7ft" 1 Mpc < a, it < 40ft~ 1 Mpc when con- 
sidering the 10k mocks and 2QZ data to minimise the risk 
of the small scale peculiar velocities affecting our ability to 
detect the geometric distortions. 

We predict that the constraints possible from £(a, it) 
measured from the 10k 2QZ data will not be as strong as for 
the finished survey, however, when we combine the results 
with the results from the amplitude of QSO clustering, we 
still manage to recover the simulation cosmology, although 
the errors are slightly larger than for the prediction from the 
finished survey. The combined results give best fit values of 
Q,^{Q)=Q.2ht° 2 2 PQSo{z)=QA0t° \l, shown in Figure 0(a). 
We just show the results assuming the A cosmology for this 
test in Figure [7| 

We also test how much the redshift errors may be af- 
fecting the constraint on cosmology. We take the positions 
of the mock QSOs from the Hubble Volume with the current 
mask applied and add a 5/i _1 Mpc Gaussian dispersion to the 
line-of-sight distances and repeat the analysis. The results 



are shown in Figure M(b) . The constraint is slightly weaker 
than when we just apply the window but is still consistent 
with the values from the Hubble Volume. The best fit val- 



-,+0.20 



We model the 



ues are fi m (0)=0.20l££; /3 QS o(z)=0.40 
peculiar velocities using an exponential model, which has 
been found to be a better fit to N-body simulations (e.g. 
Cole, Fisher & Weinberg 1994). However, if the velocities 
are dominated by redshift errors, a Gaussian model may be 
more appropriate. The fact that we are able to recover the 
Hubble Volume parameters using the exponential model even 
though we have added a Gaussian component to the small 
scale peculiar velocities suggests that the choice of model for 
the small scale peculiar velocities is not crucial for placing 
constraints on cosmology and Pqso{z). 



8.2 Results from the 10k Catalogue 

We show £(cr, 71") measured from the 10k catalogue in Figure 
pi shown by the solid lines. The dashed lines show mod- 
els, discussed below. There is evidence in the £(cr, it) plots 
that the velocity dispersion is indeed higher than 400 km 
s _1 as the data appears more elongated on small a scales 
than £(<?, it) from the Hubble Volume. We have therefore 
considered models with < w\ > 1/2 of 400, 800 and 1200 km 
S _1 . There is very little difference between the values of x 2 
found using each of the values. We obtain the best fit with 
< u>1 > 1/,2 =800 km s _1 but the constraints on cosmology 
are very similar for each value, especially when the results 
are combined with the evolution of clustering constraints. 
It is possible that the redshift errors will be reduced by the 
time the final sample is complete but for now, we only fit 
the models to the data on scales larger than 7/i -1 Mpc to 
limit the effect of the small scale peculiar velocities on the 
model. The models in Figure H and the fits shown in Figure 



are done assuming < w z > 



172 _ 



=800 km s~ 



Constraining Cosmology from Redshift Space Distortions via £(cr, n) 13 



The favoured model, if the A cosmology is assumed, 
is Q m (0)=0.3 Pqso(z)=0.35. If the EdS cosmology is as- 
sumed, then £7 m (0)=0.35 /3qso=0.45 is the favoured model. 
However, as can be seen from Figure M, the results are 
highly degenerate and a wide range of models are allowed. 
We do, however, constrain Pqso(z) to be in the range 
0.1< ,0qso < 0.5 at the la level. Currently, the only con- 
straint we can place on cosmology from fitting £(<t,7t) alone 
is that fl m (0) = is marginally rejected. All other values of 
fi m (0) are contained in the la contour. 

If we combine the results with the evolution of cluster- 
ing (Figure 0) we find the results marginally favour low val- 
ues of f2 m (0). The dashed lines show the joint la contours, 
found by adding the two values of the standard deviation to- 
gether in quadrature. If the A cosmology is assumed then we 
obtain a value of fi m (0)=0.3±g:?§, /3qso=0.35±^!?. In this 
case flm(0) — l cannot be ruled out. If we assume EdS then 
slightly lower values of Q m {0)=0.20t ;i , /3QSO=0.35±g;;[o 
are favoured. These two models are shown by the dashed 
lines in Figure tn (a, b). 

Outrarn et al. (2001) have carried out a similar anal- 
ysis of the 10k 2QZ data. They have analysed the power 
spectrum measured parallel and perpendicular to the line 
of sight and they find that the best fitting results are 
n m (0)=0.23±£^ and l3Qso=0-39t° i 8 7 - Their results are 
very similar to those found here. 



9 CONCLUSIONS 

We have developed a model for £(<r, tt) that takes into ac- 
count the different types of distortions that are introduced 
into the clustering pattern from a redshift survey, such as 
small scale peculiar velocities, bulk motions and distortions 
due to the assumed cosmology differing from the underlying 
cosmology. We have checked that the analytic models are 
consistent with £(a, tt) measured from the mock Hubble Vol- 
ume catalogues and find that the model fits the simulation 
well over the range of scales 4 < tt, a < 40/i~ 1 Mpc. 

We fit the models to £(er, tt) measured from the Hubble 
Volume and confirm previous results that cosmology, f2 m (0) 
and J7a(0), will not be strongly constrained from the 2QZ 
Survey using geometric distortions in £ (a, tt). The errors are 
such that we cannot differentiate between the effects of A 
and I3q SO (z). 

Although cosmology cannot be directly obtained from 
the measurement of £(a, tt) from the 2QZ Survey, constraints 
on /3qso(z) are possible from fitting models of t;(a, tt) to 
measurements of £(<r, tt) from the mock catalogues. This pro- 
vides a joint constraint on the cosmology and the QSO-mass 
bias. 

When these redshift space distortion results from £(u, tt) 
are combined with information on the QSO-mass bias at 
z = 1.4 from clustering evolution, new and direct constraints 
on cosmology become available. The combined constraint 
immediately rules out values of fi m (0) =0 and £7 m (0) = l with 
2a confidence, if we start with an f2 m (0)=0.3, f2A(0)=0.7 
simulation. 

We apply the methods to the 10k 2QZ catalogue and 
obtain encouraging results, although with only 10,000 QSOs 
currently any constraints on cosmology are limited. This 





Figure 9. Constraints on cosmology from £(<r, tt) measured from 
the 2QZ 10k data (greyscale) and from the evolution of the QSO 
clustering amplitude (solid lines) assuming the two different cos- 
mologies, EdS (a) and A (b). The dashed lines show the joint la 
contours. 



analysis is consistent with flat models with 0.2< il m (0) <1 
and rejects flat models with fiA(0) = l. 
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